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It is known that quantum correlations exhibited by a maximally entangled qubit pair can be 
simulated with the help of shared randomness, supplemented with additional resources, such as 
communication, post-selection or non-local boxes. For instance, in the case of projective measure- 
ments, it is possible to solve this problem with protocols using one bit of communication or making 
one use of a non-local box. We show that this problem reduces to a distributed sampling problem. 
We give a new method to obtain samples from a biased distribution, starting with shared random 
variables following a uniform distribution, and use it to build distributed sampling protocols. This 
1^-^ ' approach allows us to derive, in a simpler and unified way, many existing protocols for projective 

' measurements, and extend them to positive operator value measurements. Moreover, this approach 

, naturally leads to a local hidden variable model for Werner states. 

(N 

I. INTRODUCTION 

■ At the advent of quantum mechanics, some physicists were puzzled by the strange properties of quantum systems, 
compared to classical physics, such as randomness and non-locality. Einstein, Podolsky and Rosen showed that when 
two parties, say Alice and Bob, share an entangled state, the outcome of a measurement on Alice's side is not only 

^ ' probabilistic, but may also be conditioned on the outcome of a distant measurement on Bob's side. Therefore, they 
, questioned whether "the quantum- mechanical description of physical reality could be considered complete" 0. 
' To resolve this paradox, later called the EPR paradox, it was argued that the apparent randomness in quantum 
experiments could actually come from unknown "hidden" variables created locally along with the supposedly quantum 
state, and that this randomness would disappear as soon as these hidden variables were revealed. However, John 
Bell showed in 1964 that the quantum correlations exhibited by the EPR gedanken experiment, as reexpressed by 
Bohm f^, could not be reproduced by so-called local hidden variable models, that is, models where Alice and Bob 
share an infinite amount of locally created hidden variables Q. This proved in a sense the non-local character of 
quantum mechanics. 

Recently, in order to understand quantum non- locality, with the help of the framework of communication complexity, 
another approach has emerged to gauge the non-locality of quantum mechanics with respect to local variable models. 

■ We know that if Alice and Bob share only a set of local hidden variables (shared randomness), they cannot reproduce 
^ ' quantum correlations, but if they are allowed to use an additional resource, it may be possible for them to reproduce 

the quantum correlations. The amount of additional resources used allows us to measure non-locality. 
^ The most obvious resource that Alice and Bob can use in addition to shared randomness is classical communication. 

• I— I , The problem of simulating the quantum correlations with classical communication was studied as early as 1992 by 
' Maudlin who showed that a finite amount of communication was enough, at least on average. More precisely, 
^ , in the case of projective measurements in the real plane on a maximally entangled qubit pair, he gave a protocol 
- - ' using 1.17 bits of communication on average, but unbounded communication in the worst case. This problem was 
revisited by Brassard et al in 1999 @, who showed that 8 bits were sufficient in the worst case for arbitrary projective 
measurements. These protocols use an infinite amount of shared randomness, and indeed Massar et al. proved that 
the communication complexity can be bounded in worst case only if the amount the shared randomness is infinite. In 
2000, inspired by Feldmann 7], Steiner, independently of Maudlin, showed that for projective measurements in the 
real plane, 1.48 bits were enough on average Cerf, Gisin and Massar |^ proved that for an arbitrary projective 
measurement 1.19 bits of communication sufficed on average. Finally in 2003, Toner and Bacon lOj showed that a 
local hidden variable model supplemented with one bit of communication in the worst case is enough to reproduce 
the quantum correlations of the singlet for arbitrary projective measurements. 

Another resource which Alice and Bob can use to reproduce the quantum correlations is post-selection. Here, 
Alice and Bob are allowed to produce a special outcome of their measurement, noted _L, meaning "no result". This 
corresponds to the physical situation where Alice and Bob's detectors are partially inefficient and sometimes do 
not click. In 1999, Gisin and Gisin, inspired by Steiner's communication protocol, gave a protocol which simulates 
quantum correlations with shared randomness and with a probability 1 /3 of aborting for either party. 

Finally, Cerf et al. O have shown that a third resource could be used to simulate the quantum correlations: a 
non-local box. The non-local box is a primitive shared between Alice and Bob, with two inputs and two outputs, 
where the outputs (conditioned on the inputs) are maximally non-local in the sense that they violate a Bell inequality 
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(CHSH) maximally while remaining causal. Cerf et al. [TJ have shown that only one use of a non-local box suffices 
to simulate quantum correlations. 

In this paper, we show that the the problem can be reduced to a distributed sampling problem. In the local hidden 
variable model for two parties, Alice has an input a and an output A, and similarly. Bob has an input h and an output 
B, and they share a set of random variables which are distributed independently of Alice and Bob's input. Following an 
idea introduced by Feldmann , we can relax the condition that the shared randomness is distributed independently 
of the input, and imagine that Alice and Bob share a set of random variables with a distribution depending on Alice's 
input. Clearly in this scenario, there exists a distribution which allows them to reproduce quantum correlations (a 
trivial way is to let the random source produce a with probability 1). So the problem of reproducing the quantum 
correlations with different resources can be reduced to the problem of Alice and Bob agreeing on a sample from a 
distribution depending on Alice's input. We propose a method to carry out this distributed sampling in two steps. 
The first is the completely local problem of how Alice can sample a biased distribution depending on her input 
with the help of a (shared) uniform random source, and the second step is how Alice can share this biased sample 
with Bob with the help of various resources. After giving a new method to perform the local sampling, we will see 
that the second problem becomes easy, and allows us to understand how the various resources come into play. We 
reformulate previous protocols within this framework in an intuitive and coherent way, including the best protocols 
using communication, post-selection and non-local boxes to simulate the quantum correlations for projective or POVM 
measurements, where we also extend previous results to protocols using post-selection and non-local boxes. 

The paper will be organized as follows: in Section^ we will recall the EPR experiment and the LHV model, and 
we extend it to a setting where the parties share a biased random source. In Section [llll we will present two methods 
to perform the local sampling. In Section IIVI we will study the bipartite problem of how Alice can share her biased 
sample with Bob. First, we will see the case where Alice and Bob do not use any resource, and see that this results 
in a protocol to simulate the projective measurements on a Werner state. Second, we will study the case where Alice 
and Bob use different resources, communication, post-selection, or non-local boxes in order to simulate projective 
measurements on the singlet state. Finally, we study the previous problem for generalized measures (POVM). 



II. SIMULATING THE QUANTUM CORRELATIONS 
A. Quantum correlations 

Let us recall Bohm's version of the EPR gedanken experiment (see Fig.^: 

Definition 1 (EPR Experiment). Two parties, Alice and Bob, share a qubit pair in the singlet state, \ip^) = ( 
I Ti) ~l iT))/V2, that is, a maximally entangled state of two qubits. Alice and Bob then each receive the classical 
description of a projective measurement they have to perform on their respective qubit. These can be represented 
by unit vectors a and b pointing in some direction on the Bloch sphere. They then obtain measurement outcomes 
A e {1,-1} and B e {1,-1} respectively, where +1 corresponds to a spin parallel to the measurement direction, and 
— 1 to a spin anti-parallel to this direction. 

According to quantum mechanics, the outcome of Alice's and Bob's measurements, A and B, have the following 
joint probabilities: 

, . ^^ 1 - ^-S a • 6 
piA,B) = , (1) 

or, equivalently, their joint and marginal expectation values are given by 

E{AB\d,b) = -d-b, (2) 
E{A\d,b) = 0, (3) 
E{B\d,b) = 0. (4) 



B. Local hidden variable models 



As pointed out by Einstein, Podolsky and Rosen, the correlations between Alice's and Bob's measurement outcomes 
in such a gedanken experiment show a disturbing property of quantum mechanics, through a kind of influence at a 
distance between Alice's and Bob's outcomes, this is the EPR paradox P|. To circumvent this paradox, the following, 
completely classical, model was proposed to simulate the EPR experiment (see Fig.|2l 



Alice 
A 



Bob 
B 



FIG. 1: EPR experiment. Alice and Bob share a pair of qubits in the singlet state lip ) = (| TJ.) ~ 1 11))/V2. Both perform 
measurement on their qubit, specified by vectors a or b, and obtain results A — ±1 or B — ±1. 



Definition 2 (Local hidden variable model). Alice and Bob share some random variable A G A, where A is some 
possibly infinite set, with probability distribution p{X). They then receive inputs a G §2 and 6 G S2 respectively, and 
output A — A(d, A) G {1, —1} and B = B{b, A) G {1, —1} respectively. 

In this scenario the functions A{a, A) and B{b, A) are deterministic so that, as intended, Alice's and Bob's outputs 
are fixed as soon as the value of the random variable A is known. Moreover, this hidden variable may have been 
created locally (as was the quantum state in the EPR experiment) and then communicated to both Alice and Bob. 
Therefore, if such a model could reproduce the quantum correlations of the EPR experiment described above, it could 
conveniently replace quantum mechanics as it would circumvent not only randomness but also non-locality. 
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FIG. 2: LHV model for the EPR experiment. 

However, Bell has shown in his famous theorem that such a simulation is not possible 0. Indeed, let us note that 
we have implicitly assumed that the hidden variable A was distributed independently of the measurement directions a 
and 6, p{X\d, b) = p(A), since the physical interpretation is that A has been created along with the supposedly quantum 
state, which could have happened long before the inputs a and b were fixed. 

Theorem 3 (Bell). No local hidden variable model may simulate the quantum correlations exhibited by the EPR 
experiment given in Definition^ 

Proof. (Sketch) For a local hidden variable model, we have E{AB\a,b) = T^gaPW A{a, X)B{b, X). It is then easy 
to show that all local hidden variable models satisfy the CHSH inequality 



where 



^2 < C(ai, a^2, bi,b2) < 2 (Vol, 61, 0^2, &2 e §2), 



C(ai, a^2, bi,b2) = E{AB\ai,bi) + E{AB\a-{,b2) + -B(ylS|al, fe^) - £;(AB|al, 62) 



(5) 



(6) 



Nevertheless, for quantum mechanics, there exist dl, 02, bi, &2 G §2 such that Eq. ^ yields C{di, 02, bi, 62) ~ 2^/2 > 2, 
so that quantum mechanics violates the CHSH inequality and therefore cannot be reproduced by a local hidden variable 
model. □ 



C. Protocol with a biased random source 



We have seen that a local hidden variable model, as defined above, where the random variable A was distributed 
independently of the inputs a and b, does not allow for the simulation of the EPR experiment. Nonetheless, if we relax 
our model and let the random variable depend on one of the inputs, it becomes possible to reproduce the quantum 
joint distribution ^ (this is a slight extension of a result by Feldmann 



4 



Theorem 4 (Sampling theorem). Let a and b be Alice's and Bob's inputs. If Alice and Bob share a random 
variable Xg G S2 distributed according to a biased distribution with probability density 

••a: 

p{Xs\S,b) ^ PsiXs) ^ — , (7) 



2tt 

then they are able to simulate the EPR experiment in Definition^ without any further resource, that is, simulating 
the EPR experiment reduces to distributed sampling from p^. 

Proof. If Alice and Bob set their respective outputs as A(a, As) = — sgn(a • As) and B{b,Xs) = sgn(6 • As), where 
sgn(x) = 1 for a; > and sgn(a;) = —1 for x < {x G K), then the joint expectation value E{AB\a, b) is given by 

E{AB\d,h) ^ [ p{Xs\d,b) A{d,Xs)B{b,K) dK (8) 

JS2 

= -IT [ {a ■ K) sgnib ■ X,) dXs (9) 
Js2 

= -a-b. (10) 

Similarly, we have E{A\d, b) = E{B\d, b) = 0, as desired. □ 

Theorem 0] allows us to reduce the simulation of quantum correlations to the distributed sampling of a biased 
distribution p^. The problem can further be reduced to two steps. First, Alice can locally create a sample As ~ ps 
using her knowledge of a and uniformly distributed random variables. This is what we will call local sampling. The 
second step is a communication complexity problem: Alice shares this biased variable with Bob. For this step, they 
will need resources in addition to those allowed by local hidden variable models. Let us now study the first problem, 
the local sampling. 



III. LOCAL SAMPLING OF THE BIASED DISTRIBUTION 



A. The rejection method 

From now on, we will use the notation As for the biased samples, (Aq, Ai, . . .) for the sequence of uniformly distributed 
random variables shared by Alice and Bob, and Ua for a uniform distribution on the set A. 

Forgetting about Bob, the problem is for Alice to sample from the biased distribution pg from a source of uniformly 
distributed variables. A well known method to perform (local) sampling is the rejection method [l3j |. which in our 
case gives: 

Theorem 5 (The rejection method). Let Alice perform the following protocol: 







Rejection method 


Set k 


= 




1. 


Alice picks Xk 




2. 


Alice picks Uk 




3. 


Ifuk < \a- Xk\ 


, then she accepts Xk and sets Xs — Xk, 




otherwise, she 


rejects Xk and goes back to step 1 with k = k + 1. 



When the process terminates, we have As ~ |a • As|/27r. 

Let us note that this method is iterative: if Alice rejects A^, she has to start over with a fresh random variable 
Afc+i, and so on until she accepts one sample. If Alice is particularly unlucky, she could reject an arbitrarily large 
number of samples before accepting one. This means that in the worst case the process takes unbounded time to 
terminate and requires an infinite amount of samples. 
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B. The "choice" method 

To avoid this drawback, we now propose a new method, that we wiU call the "choice" method. Contrary to the 
rejection method, this method will not reject anything, but it is less general than the rejection method because it 
makes use of a specific property of the biased distribution that we want to sample from. Indeed, the key is the 
following remark: in the rejection method above, the bias \a ■ Xq\ is uniformly distributed in [0,1] when Xq ^ Ug^. 
Therefore, we could produce a sample uq ^ ^[o,i]: by picking a second vector uniformly distributed on the sphere 
Ai Ugj, and setting uq = \d ■ Ai|. Hence, to sample from the biased distribution ps{Xs) = |a ■ As|/27r, we may use 
the following theorem: 

Theorem 6 (The "choice" method). If Alice performs the following protocol: 

Choice method 

1. Alice picks Aq ^ Ug^, 

2. Alice picks Ai ~ Ug^ , 

3. //|a- All < |a- Aol, then she accepts Ao and sets Xg — Xf), 
otherwise, she accepts Ai and sets Xg = Ai. 

then Xg ^ \a, ■ Xs\/2tt and p{Xs = Aq) = p{Xs = Ai) = ^. 

Proof. Let uq — |a • Ai| and ui ^ \a ■ Ao|. Then mo,ui ~ '-'[0,1] when Aq, Ai ^ Ug^, so the rejection method ensures 
that conditionally to the test uq < |a- Ao| being satisfied, A^ will be distributed according to the distribution p^. 
Exchanging Aq and Ai's roles, the same argument ensures that conditionally to the test ui < |a • Ai| being satisfied. 
As will also be distributed according to ps, so we have p{Xs\uo < \a ■ Ao|) — p{Xs\ui < \a ■ Ai|) = ps{Xs). As the 
outcomes of the test are of course mutually exclusive. As will always have the right distribution. □ 

IV. DISTRIBUTED SAMPLING AND SIMULATING CORRELATIONS 

A. Simulation of a Werner state 

A Werner state is a quantum state obtained by mixing a pure state \ip) with maximally random noise. It is 
characterized by a density matrix W — p\ip){ip\ + (1 — p)ll/(i, where d is the dimension of the Hilbert space in which 
\ip) lies, and p is the so-called visibility of {ip) in the Werner state. These states were introduced by Werner to show 
that the correlations exhibited by some non-separable states could be simulated by local hidden variable models [T^ . 
We will now show that our approach leads naturally to a local hidden variable model for a Werner state of the singlet 

IV'->. 

In the previous section, we have given two methods for Alice to sample from the specific biased distribution using 
uniformly distributed variables. Suppose Alice performs the choice method. Without any further resource than that 
allowed by local hidden variable models. Bob is not aware of Alice's choice, and hence does not share with Alice the 
variable As ^ ps- Because of Bell's theorem, we know that Alice and Bob will not be able to simulate the singlet 
correlations, but Bob could nevertheless try to guess As. More specifically, if he assumes that As = Aq, he will be 
right half of the time. As an intermediate step towards simulating singlet correlations, we show that with such a local 
hidden variable model, Alice and Bob may simulate a Werner state of the singlet state with a visibility of one half. 

Theorem 7 (Simulation of Werner states). A local hidden variable model may simulate the correlations exhibited 
by the Werner state 

W^p\i;-){4'-\ + {l-p) ^ (11) 

with a visibility of p = 1/2. 

Proof. As already stated, the local hidden variable model consists for Alice to perform the choice method, and for 
Bob to always assume that Alice chose Aq: 
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LHV model for Werner state 

1. Alice and Bob share a pair of uniformly distributed variables Aq, Ai ~ Ug^, 

2. Alice performs the choice method: she tests whether |a • Ai| < |a • AqI, 

• if so, she outputs A = — sgn(a • Aq), 

• otherwise, she outputs A = — sgn(a • Ai), 

3. Bob outputs B = sgn(6 • Aq). 



To prove that this model achieves its goal, we first note that Alice's output may be rewritten in the more compact 
form A = — sgn(a • (Aq + Ai)). The local hidden variable then yields: 

E{AB\d,b) = I dXo I dA^i sgn(a- (A'^ + AO) sgn(6- A"o) (12) 



(47r)2 
a ■ b 



(13) 

that is the same joint expectation value as for the Werner state (lllll with a visibility of p = 1/2. As we also have 
E{A\a, h) = E{B\a, h) = 0, this model reproduces all the correlations exhibited by this state. □ 

Actually, even though we have derived our local hidden variable in a quite different way, it finally reduces to a 
rewriting of the model proposed by Werner in p^ . 

B. Simulation of the singlet state with additional resources between Alice and Bob 

Now, if Alice and Bob want to simulate the correlations of the singlet with perfect visibility, we know from Bell's 
theorem (Theorem (S)), that they will need additional resources. Moreover, we know that they may achieve their 
goal as soon as they share some random variable As distributed according to the biased distribution p^. We will 
now consider three different resources that will allow Alice and Bob to sample from this distribution: post-selection, 
classical communication, and non-local boxes. Let us note that in this section, we will not introduce any new protocols, 
but rather we will show that the most efficient known protocols can be viewed as sampling from the same biased 
distribution ps- This reduction to the distributed sampling problem is our main contribution. 



1. Post-selection 



Let us first consider post-selection, which is the possibility for Alice and Bob to sometimes abort the protocol. 
Physically this would correspond to the fact that Alice and Bob's detectors are imperfect and sometimes do not click. 
In our model, post-selection is represented by a new symbol _L added to the set of valid outputs for Alice and Bob. 
We will use the notation J- a for the event that Alice aborts, and similarly _Lb for Bob. 

Theorem 8 (Post-selection). If Alice and Bob are given an infinite amount of shared randomness, supplemented 
with post-selection for only one of them, say Alice, there exists a protocol with p(_L^) = 1/2 such that conditionally 
to the fact that no party aborted, Alice and Bob share a random variable A^ ~ ps. 

Together with Theorem 0] this implies the following corollary 

Corollary 9. Using a local hidden variable model supplemented with post- selection, it is possible for Alice and Bob 
to simulate the EPR experiment in Definition^with p{J-a) = 1/2 and p{J-b) — 0. 

Proof of Theorem\^ The idea is to let Alice perform the rejection method (Theorem (S)) and use post-selection to 
reject the bad A's, which leads to the following protocol: 
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Post-selection protocol 


1 

J_ . 


r\liCe dliLl JDOU bildl c d UlillOl iliij CllbLilUULeU. VdiiclUic A 1 ^§2 1 


z. 


Alice picks u ^ L-l[o^i], 


3. 


She perforins the rejection method: she tests whether u < \d ■ A|, 




• if so, she sets As — A, 




• otherwise, she aborts the protocol {-La), 


4. 


Bob always sets As — A. 



Since both u and |a • A| arc uniformly distributed on [0, 1] when A ~ Ug^ , it is clear that p{-La) — p{u < |a • A|) ~ 1/2. 
Moreover, Theorem |S1 ensures that p{Xs\^±a) — Psi^s)- □ 

Let us note that this is just a rewriting of Gisin and Gisin's protocol Here, only Alice uses post-selection, 

but, as pointed out by these authors, it is clearly possible to symmetrize the protocol by randomly exchanging Alice's 
and Bob's role, and using an additional shared random variable to tell Alice and Bob to both abort simultaneously 
(see ^3 fo'" details). This leads to a protocol with p{-La) — p[^b) = 1/3. 

2. Communication 

Let us now consider a second resource used to sample from p^: communication between Alice and Bob. First of all, 
we briefly recall Steiner's protocol [sl, which, in a similar manner to Gisin and Gisin's, consists in using the rejection 
method to sample from the biased distribution. His idea was to consider that Alice and Bob share an infinite sequence 
(Ao, Ai, • • • ) of samples A^ ~ Ug^ (recall that the local hidden variable model assumes no limitation on the amount 
of shared randomness). Similarly to Gisin and Gisin's protocol, Alice then performs the rejection method with the 
first sample Aq, but instead of aborting the protocol if she has to reject Aq, she iterates the method by taking the 
samples in the order of the sequence until she accepts one of them, say A^ . She then communicates the index k of the 
accepted sample to Bob. As pointed out earlier, if Alice is particularly unlucky, she could reject an arbitrarily large 
number of A's before accepting one. The index k she has to communicate to Bob may then become arbitrarily large, 
such that in the worst case the amount of communication required to simulate the EPR experiment with this method 
is unbounded. 

Instead of using the rejection method, which leads to unbounded communication complexity in the worst case, we 
can use the choice method to obtain a protocol with bounded communication in the worst case. 

Theorem 10 (Communication). If Alice and Bob are given an infinite amount of shared randomness, supplemented 
with one-way communication, there exists a protocol using exactly one bit of communication that allows Alice and Bob 
to share a random variable Xg ps- 

Together with Theorem 0] we have the following corollary: 

Corollary 11. Using a local hidden variable model supplemented with one-way communication, it is possible for Alice 
and Bob to simulate the EPR experiment in Definition^ with exactly one bit of communication. 

Proof of Theorem MIA Let the two parties use the following protocol, where Alice performs the choice method shown 
in Theorem El and then communicates to Bob the index of the sample she accepted: 
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Communication protocol 


i. 


Alice and Bob share a pair of miifornily distributed variables Aq, Ai ^ Llgj, 


2. 


Alice performs the choice method: she tests whether a - Ai < a • Aq , 




• if so, she accepts Aq and sets a; = 0, 




• otherwise, she accepts Ai and sets x = 1, 


3. 


Alice sends x to Bob, 


4. 


Alice and Bob then set A^ = Ax- 



Hence, Theorem ensures that As ~ |a • Xs\/2tt. □ 

Let us note that CoroUarv lllI was proven by Toner and Bacon p^ . but our approach has the advantage of giving 
much more intuition, and in particular of clarifying the relationship with the previous protocols. 

3. Non-local boxes 

We now consider another additional resource shared by Alice and Bob: a non-local box. This resource was introduced 
by Popescu and Rohrlich , and has recently been used by Cerf, Gisin, Massar and Popescu to simulate the singlet 
correlations Let us recall its definition and its main characteristics. 

Deflnition 12. A PR non-local box is a device shared by Alice and Bob, that has two input bits x,y € {0, 1} from Alice 
and Bob, respectively, and outputs a, P € {0, 1} to Alice and Bob, respectively, according to the following distribution: 

I M \ \ \ ifa®P = xAy, , . 

p(a,/3|x,y) = | 2 ^^^^^^^^^^ (14) 

One use of a non-local box will be called an nl-bit. This resource has the following interesting properties: 

• it is maximally non-local, in the sense that it maximally violates the CHSH Bell inequality. 

• it is causal, in the sense that Alice's output a is independent of Bob's input y, p{a\x,y) — p{a\x) (and vice 
versa) , 

• it is a strictly weaker resource than one bit of communication: due to the causality property, it may not be used 
to communicate but, on the other hand, it may be shown that one use of a non-local box may be simulated by 
one bit of communication jil| . 

Since we have shown that the previous protocols for simulating the singlet correlations reduced to sampling from 
the distribution pg, and motivated by Cerf, Gisin, Massar and Popescu's result 11], we could be tempted to study 
whether a non-local box also allows Alice and Bob to share a random variable A^ ~ p^. However, this is obviously 
not possible since As is not independent from a, so that their mutual information is non-zero, /(As : a) ^ 0. Indeed, if 
Bob knew As, this would mean that he has gained some information about a, but this is impossible using only shared 
randomness and a non-local box since we have seen that a non-local box does not allow signaling. 

However, by examining the proof of Theorem^ we notice that to simulate the singlet correlations from a shared 
random variable As psi Bob actually only uses the knowledge of sgn(6 -As), (and similarly for Alice, but this is not 
really an issue). Let us introduce the following definition: 

Definition 13. Let d,b be Alice's and Bob's inputs. A protocol is an f -sampling protocol for a given distribution p 
if at the end, Alice has a sample sa = /(a. As) and Bob has a sample sb = f{b,Xs), where As ^ p. 

We may now rewrite Theorem0]with a weaker sampling hypothesis: 

Tiieorem 14 (Strong sampling theorem). The simulation of the EPR experiment in Definition^is eguivalent 
to f -sampling for pg, where ps(As) = |a- As|/27r, and f(x,y) = sgn(x ■ y). 
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Proof. The proof that simulating the EPR experiment reduces to /-samphng for pg is identical to Theorem 0] To 
show that /-sampling for reduces to simulating the EPR experiment, let A, B be the outcome of a simulation of 
the EPR experiment. Let sa = —A and sb = B. Then sa, sb is an /-sample for p^. 

□ 

With this weaker hypothesis in mind, we will now prove the following theorem: 

Theorem 15 (Non-local box). If Alice and Bob are given an infinite amount of shared randomness, there exists 
an f -sampling protocol for ps with f{x, y) — sgn{x ■ y) that only makes use of 1 nl-bit. 

Once again, together with Theorem 1141 we have the following corollary: 

Corollary 16. Using a local hidden variable model supplemented with a shared PR non-local box, it is possible for 
Alice and Bob to simulate the EPR experiment in Definition^ with exactly 1 nl-bit. 

Proof of Theorem \15\ Before introducing the protocol, let us give some motivation by recalling the protocol that 
simulates the Werner state with visibility p — 1/2 (see Theorem II If) . In this protocol, Alice performs the choice 
method (she chooses between Aq and Ai) and Bob always assumes that she chose Aq. Obviously, he will be wrong half 
of the times, when Alice chose Ai, but sometimes this is not a problem since it may happen that sgn(5-Ao) — sgn(6- Ai). 
The idea of the following protocol is to use a PR non-local box to correct Bob's remaining mistakes, which happens 
when Alice chooses Ai and sgn(6 • Aq) 7^ sgn(6 • Ai). 





Non-local box protocol 


1. 


Alice and Bob share a pair of uniformly distributed variables Aq, Ai ~ Ug^, 


2. 


Alice performs the choice method: she tests whether |a • Ai| < |a • Ao|, 




• if so, she accepts Aq and sets a; = 0, 




• otherwise, she accepts Ai and sets a; = 1, 


3. 


Alice inputs x into the non-local box. 


4. 


Bob tests whether sgn(6, Aq) = sgn(6, Ai) 




• if so, he sets y = 0, 




• otherwise, he sets y — \, 


5. 


Bob inputs bit y into the non-local box. 


6. 


Alice gets output bit a from the box and sets As — (— 1)"A3;, 


7. 


Bob gets output bit (3 from the box and sets A^ = (— l)'^Ao. 



Let us now analyse the protocol. As Alice uses the choice method to choose Xx, we know from Theorem El that 
Xx Pa- Moreover, since a is an unbiased random bit, As = (— l)"Aa; is also distributed according to p^. It is clear 
that Alice may compute sa — sgn(a • As) since she knows a and Ag. 

On Bob's side, an error takes place if and only if Xx = Ai and sgn(fe, Aq) 7^ sgn(5, Ai), and this is precisely the case 
when the non-local box's outputs are different: 

[Xx = Ai and sgn(6, Aq) ^ sgn(6, Ai)] -^=» [x A y] -^=^ /3 ^ a. (15) 

Therefore, the non-local box corrects Bob's error by changing the sign of Xb = (— l)''Ao with respect to As when 
necessary, so that he may correctly compute sb — sgn(6 ■ Xb) = sgn(5 -As). □ 

Let us note that this reproves Cerf, Gisin, Massar and Popescu's result pTl |. but with a more intuitive approach. 



10 



C. Simulation of POVMs with additional resources 



In the last section, we derived protocols for simulating local projective measurements on the singlet. However, 
quantum mechanics allows a wider class of measurements than projective measurements, and we will now use our 
approach to simulate the EPR experiment on the singlet with the most general type of measurement allowed by 
quantum mechanics, namely positive operator value measure (POVM). Let us first recall some definitions. 

Definition 17 (POVM). A positive operator value measure A of size r on a Hilbert space 7i„ of dimension n is a 
collection of r positive (i.e. with non-negative eigenvalues) operators on 7i„, {Ai\i = 1, . . . ,r}, such that 

r 

J2a,^1L, (16) 

i=0 

where lln is the identity on 7i„ . 

The projective measurements that we have considered above are a particular case of POVM where the elements Ai 
are orthogonal rank one projectors, that is AiAj — SijAi. In other words, they define a basis {li^i)!? = 1, . . . ,r} such 
that A, = 

Without loss of generality, we can restrict our study to rank one POVMs, i.e., POVMs where the elements Ai are 
proportional to rank one projectors \(f)i){(j)i\. In dimension 2, that is for a qubit, a rank one projector \4>i){(t>i\ may be 
represented on the Bloch sphere by a unit vector v such that \(l)i){4>i\ = {II2 + v ■ (t)/2, where a is the vector of Pauli 
matrices. Therefore, the elements Ai of a rank one POVM may be represented by a (not necessarily unit) vector a; 
such that Ai — (|ai|ll2 + Si ■ (j)/2. 

From Definitional we see that a rank one POVM A of size r on a qubit may be represented by a set of r vectors 
in the Bloch ball, {di\i = 1, . . . ,r} satisfying the following conditions, equivalent to 

r 

El«^l = 2, (17) 

r 

^al = 0. (18) 

1=0 

Now we generalize the EPR experiment in Definition to the case where Alice and Bob perform POVMs. 

Definition 18 (EPR Experiment with POVMs). Two parties, Alice and Bob, share a qubit pair in the singlet 
state, — (I tl) — I |'f))/\/2. Alice and Bob each receive the classical description of a POVM they have to 

perform on their respective qubit. These can be represented by collections of vectors A = {di\i — 1, . . . , rA} and B = 
{bj\i = 1, • . • , fs} satisfying conditions jiTI j<*i|] . They then obtain measurement outcomes Apov G {1, r^} 
and Bpov G {1, . . . , j, . . . , rs} respectively. 

According to quantum mechanics, the outcomes of Alice's and Bob's measurements, Apov and -Bpov, have the 
following joint probabilities: 

p(Apov = z,i?pov=jiAs) = '°'"^^^'7"''•^^ (19) 



with marginal probabilities 



p(Apov = *|AS) = ^, (20) 



2 ' 

2 



p(Spov=j|-4,S) = M. (21) 

(22) 



We wish to simulate these correlations between Apov and Bpov using a local hidden variable model with the help 
of additional resources. We have seen above that it is possible to simulate projective measurements on the singlet 
by sampling or /-sampling a random vector according to a biased distribution psi^s) — \a ■ Xs\/2tt (Theorems01 
and ll4|) . that is, if Alice and Bob have as input unit vectors a'^ = di/\ci^i\ and b'j = bj/\bj\, and an infinite amount of 
shared randomness, they are able to produce, with additional resources, random variables A = ±1 and B = ±1 with a 
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joint distribution p{A, 5|a^, bj) = (1 — AB a\ ■ &j)/4. In order to simulate the POVMs, we will use one of the previous 
protocols to simulate projective measures, then we will test whether the outcomes of the projective measure protocol 
"agree" with the POVM outcome we expected. We call this test the POVM test. In addition to the resources used 
for the projective measure protocol, the POVM test will require further resources: post-selection or communication. 
In this approach, the test cannot be performed using only non-local boxes, because of their non-signaling property. 

Theorem 19. It is possible for Alice and Bob to simulate the EPR experiment for POVMs using any one of the 
following ressources. 

(a) postselection, with p{J-j^) =p(_Lb) ~ 2/3; 

(b) 6 bits of communication on average; 

(c) 2 nl-bits plus 4 bits of communication, on average. 
Proof. Consider the following protocol. 

POVM protocol 

Set fc = 

1. Alice and Bob share a pair of uniformly distributed variables A^, Afc+i ~ Ug2, 

2. Alice picks Si G A following the marginal probabilities p{i) — \di\/2. 

3. Bob picks bj e B following the marginal probabilities p{j) = \bj\/2. 

4. Alice and Bob perform a protocol to simulate the projective measurements on 
a singlet with additional resources (either post-selection, communication or a 

non-local-box) using a'^ — di/\ci^i\ and b'j — bj/\bj\ as inputs, and obtaining A 
and B as outputs, 

5. POVM test: 

• with post-selection: 

(a) If A = 1, Alice outputs Apov — i, otherwise she aborts, 

(b) li B = 1, Bob outputs i?pov — j, otherwise he aborts. 

• with communication: 

(a) Ahce sends A to Bob, 

(b) Bob sends B to Alice. 

(c) Alice and Bob test whether A — B, 

— if so, they output Apov — i, and i?pov = J respectively, 

— otherwise, they return to step 1 with k ~ k + 2. 

Intuitively, Alice and Bob choose POVM outcomes dl G A and bj € B with the right marginal distribution, and 
they use the projective measurement protocol to realize this outcome. To do this, they simulate the corresponding 
(normalized) projective measurements a'^ and b'j on the singlet. Alice (and similarly for Bob) then gets an output 

A{a'^, Xg) — ±1, corresponding to a spin either parallel or antiparallel to dl. While an outcome -1-1 corresponds to 
a valid POVM element dl e A, —1 does not correspond to a POVM element, —di ^ A, so a priori, Alice should 
only accept -1-1 outcomes. In the POVM test, they test whether the outcomes of the corresponding projective 
measurements represent valid POVM elements a'^ and b'j (POVM test), and if not, either start over again or abort if 
they use post-selection. 

Finally, depending on the resource used by Alice and Bob to perform the projective measurement simulation, we 
have the following results. 

Post-selection If Alice and Bob use post-selection, using Theorem [S] they may simulate the projective measure on 
a singlet with a symmetric and apparently independent abortion probability of 1/3 for both parties. Moreover, both 
Alice and Bob will abort when A or B equals —1, which happens with probability 1/2, so that they output only when 
the POVM test is satisfied. Altogether, they finally use post selection with probability 2/3, p{-La) — p(-Ls) — 2/3. 
This result is new. 
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Communication Recall that an outcome +1 of the projective measurement corresponds to a valid POVM element 
di & A However, considering both Alice and Bob's outcomes together, we see that A — B = —1 happen with the 
same probability that A ^ B = +1, so that in the end, Alice and Bob only have to reject their outcomes when 
A ^ B. If Alice and Bob use communication, from Theorem I lUI we see that they must use one bit of communication 
to simulate the projective measure on the singlet state, and two bits of communication to perform the POVM test. 
As the POVM test has an average probability of being satisfied of 1/2, they must repeat the projective measurement 
protocol twice on average to satisfy the test. So, they use six bits of expected communication. This is a rewriting of 
Methot's result 

Non-local boxes and communication Considering the last protocol, we may wonder whether the six bits of 
expected communication may be replaced by uses of a non-local box. Obviously, we may simulate the projective 
measure on the singlet with 1 nl-bit. However, to perform the POVM test in our protocol, two-way communication 
is necessary since to know whether her outcome A agrees with Bob's, Alice needs to know B and therefore acquire 
information from Bob, and vice versa. Hence, our POVM test cannot be done with non-local boxes only since they 
do not allow signaling. In this scenario, the best we can achieve is a protocol that uses two nl-bits and four bits of 
communication on average. This result is new. □ 



V. CONCLUSION 



In summary, we have shown that the problem of simulating quantum correlations exhibited by the singlet state 
using different resources reduces to a distributed sampling problem, more precisely, to an /-sampling problem. We 
have seen that to perform distributed sampling from a biased distribution depending only on Alice's input, Alice 
can first sample the biased distribution locally and then use an additional resource to share the biased distribution 
with Bob. To locally sample from the biased distribution pg, we have given a new method called the "choice method" . 
Once the local sampling is done, Alice uses an additional resource to share her sample with Bob in a simple manner. 
This approach allows us to develop a simple and unified view of the problem and therefore a better understanding of 
the role of the different resources used to gauge non-locality. The distribution has a crucial role in our approach 
and we may ask whether it is the unique distribution that arises naturally from this problem. We have seen that 
simulating the EPR experiment is equivalent to /-sampling from ps, but is not equivalent to distributed sampling 
from ps- An intermediate sampling problem is for Alice to sample from p^ and Bob to learn only sgn(& ■ Xg). It would 
be interesting to know whether this problem is equivalent to simulating the EPR experiment or not. Such a result 
would help us understand other properties of this problem, such as the optimality of the protocols. The POVM case 
also fits into our approach, but the protocol we have derived is probably not optimal because it simply consists in 
adapting the projective measurement protocol to the POVM case. It might be the case that protocols can be based 
on the /-sampling of another biased distribution designed especially for POVMs. 

Finally, this approach could also help study more general cases of the problem, for instance for higher dimensional 
states (such as maximally entangled qudit pairs), multipartite entangled states, or partially entangled states. 
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